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Abstract 

We study relativistic effects in a system of two scalar particles interacting via a scalar 
exchange in the Light Front Dynamcis framework. The results are compared to those 
provided by Bethe-Salpeter and non relativistic equations. It is found in particular 
that for massive exchange, the relativistic description is of crucial importance even 
in the limit of zero binding energy. 
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1 Introduction 



Light Front Dynamics (LFD) is a field theoretically inspired hamiltonian ap- 
proach specially well adapted for describing relativistic composite systems. 
First suggested by Dirac [1] it has been since widely developed (see [2-6] and 
references therein) and recently applied with success in its explicitly covari- 
ant version [7] to the high momentum processes measured in TJNAF [5,8]. In 
this approach the state vector is defined on a space-time hyperplane given by 
uj ■ x = a where uj = (1,77.) is a light-like four- vector. 

We present here the first results obtained within this approach for the Wick- 
Cutkosky model [9] . This model describes the dynamics of two identical scalar 
particles of mass m interacting by the exchange of a massless scalar particle. 
This first step towards more realistic systems constitutes an instructive case 
and is presently considered by several authors [10-14]. The model has been 
extended to the case where the exchanged particle has non zero mass \i and 
used to build a relativistic scalar model for deuteron. 



The results presented in this paper concern the S-wave bound states in the 
ladder approximation. They are aimed to (i) compare the LFD and Bethe- 
Salpeter descriptions and study their non relativistic limits, (ii) disentangle 
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the origin of the different relativistic effects, (in) evaluate the contribution 
of higher Fock components and (iv) apply this study to a scalar model for 
deuteron. 



2 Equation for Wick-Cutkosky model 

We have considered the following lagrangian density: 

C = X - (d v (j&t - m 2 2 ) + \ (duXd u X - A*Y) - 9<P\ 

where and x are re& l fields. In the case /j, — it corresponds to the Wick- 
Cutkosky model. The wave function describing a bound state of two par- 
ticles with momenta k\ and k 2 , satisfies in the Light-Front the dynamical 
equation [7,5] 

[4(q 2 + m 2 ) - M 2 Mq, n) = J ^V(q, q', n, M 2 )^{q\ h) (I) 

Variable q is the momentum of one of the particles in the reference system 
where k\ + k 2 = 0, and tends in the non relativitic limit to the usual center 
of mass momentum. M represents the total mass of the composite system, 
B = 2m — M denotes its binding energy and e q = y/q 2 + m 2 . In the case of 
S-waves the wavefunction is a scalar quantity depending only on scalars q and 
n- q [7,5]. 

The interaction kernel V calculated in the ladder approximation is given by 

W,?>,M>)=-^ (2) 

with 

Q 2 = (q- q>Y -(n.q)(n. ^^f- + (f 9 + 4 - ^) 

The coupling parameter a is defined by a = g 2 /16nm 2 . In the limit « 2m 
some analytical solutions are known [15] and once removed the n dependence 
in the kernel - formally setting n — - the LFD equation turns back to the 
Schrodinger equation in the momentum space for the Yukawa or Coulomb 
potential, a being the usual fine structure constant. 



n ■ q' n ■ q 
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Equation (1) has been solved with the coordinate choice displayed in figure 1. 
We have chosen z axis along h and with no loss of generality y? = 0. The ip' 
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Fig. 1. Kinematical variables. The z axis is chosen along n and q can be restricted 
to the xOz plane without loss of generality 

dependence of the kernel (2) can be performed analytically and (1) turns into 
the two dimensional integral equation 

[A(q 2 + m 2 ) -M 2 ]^(q,9) = f ^dq'sm9'd9' 1 V(q', 9') (3) 



with 



a = q 2 + q' 2 - qq' 2 cos 9 cos 9' + 



- £ q y 



+ [q 2 + q' 2 + 2m 2 - 



b = 2qq'sin9sm9' 



M 2 



q' cos 9' q cos 9 



The kernel of (3) has an integrable singularity for (q, 9) = (q',9 r ). The equa- 
tion is solved by expanding the solution on a spline functions basis Si, asso- 
ciated with coordinates q and 9: ^(q,9) = J2i,j CijSi(q)Sj(9). The r.h.s. two- 
dimensional integral is evaluated using Gauss quadrature method adapted to 
treat the singularity. The unknowns of the problem are the coefficients c^-, 
which are solutions of a generalized eigeinvalue problem XBC = A(M 2 )C for 
M 2 values such that A(M 2 ) = 1. 



3 Results 



The LFD binding energy for /i = versus the coupling constant is displayed in 
figure 2 (solid line). It is compared with the non relativistic values (dot-dashed 
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line) and a first order perturbative calculation B pert (dashed line), valid also 
for Bethe-Salpeter (BS) equation [16], given by 



Bpert 



ma 2 



(l + ^aloga) (4) 



Corresponding numerical values - ink — c — m — 1 units - are given in table 
1 together with the quantity R = <<? 2 > usually used to evaluate the relativistic 
character of a system. A first sight at this figure shows a significant departure 
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Fig. 2. Binding energy (solid) for fj, = compared to non relativistic (dot-dashed) 
and perturbative (dashed) calculations 

from the non relativistic results already for a — 0.1. This discrepancy - which 
keeps increasing till B reaches the maximum value of 2m - is of 100% for 
a = 0.3 whereas R remains very small. When evaluated using non relativistic 
solutions, R is equal to B (virial theorem), what gives R ~ 2% for a = 0.3, in 
contrast with the 100% effect in the binding. The R values obtained using the 
LFD solutions are even smaller (see table 1). It is worth noticing the sizeable 
relativistic effects observed in a system for which both the binding energy and 
the average momentum are small. 

A good agreement with the perturbative calculation is found up to values 
a 0.3 where the relative differences are 3%. The particular form of equation 
(4) ensures the existence of a non relativistic limit, the same for LFD and 
BS approaches, for the weakly bound states. We will see later on, that this 
situation is particular to the case fj, — 0. 



The bound-state wavefunctions presented below are normalized according to 



Table 1 

Binding energies for \i = and fj, = 0.15 as function of the coupling constant a 





0.3 


0.4 


0.5 


1.0 2.0 


3.0 


4.0 


5.0 


6.0 


6.98 


7.26 


10 2 B^q 


1.17 


1.87 


2.67 


7.68 21.0 


38.0 


58.5 


84.0 


118 


200 




10 2 ^> 


1.1 


1.8 


2.5 


6.7 17 


28 


40 


55 


64 


80 




10 2 5^=0.15 




0.190 


0.570 


4.36 16.5 


32.7 




75.6 


107 


157 


200 


with e q = m 


for the 


non 


relativistic case. The LFD 


wave 


function ^f(q,9 = 


0) 



obtained for a = 0.5 is compared in figure 3a (solid line) with the correspond- 
ing non relativistic solution (dot-dashed line), that is Coulomb wave function. 
The sizeable difference between both functions is mainly due the differences 
in their binding energies: Blfd = 0.0267 whereas B^r = 0.0625 for the same 
coupling constant. In order to compare wave functions with the same energy, 
the value of the coupling constant for the non relativistic solution is adjusted 
to ajyji = 0.327. The wave function obtained (long-dashed curve) is then much 
closer to the relativistic one. 
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Fig. 3. (a) S-wave LFD solution for 6 = (solid) compared with Coulomb wave func- 
tion with the same coupling constant (dot-dashed) and the same energy (dashed). 
Figure (b) shows the angular dependence of the wave function for different values 
of q. Curves from bottom to top correspond to increasing momenta from q = to 
q = 1.5. 

Furthermore, in the region of high momentum transfer, the relativistic func- 
tion is smaller than the Coulomb one, as expected from the natural cut-off of 
high momentum components introduced by relativity. However, these differ- 
ences can be accounted with the 6 dependence of the LFD solutions, which 
exists even for S-waves. This angular dependence, normalized by the value of 
^(q, 6 — 0), is shown in figure 3b for different values of momentum < q < 1.5. 
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As one can see, the influence of the momentum orientation compared to the 
light-front plane is far from being negligible. This effect increases with q and, 
for a fixed value of the momentum, is maximum when q ■ n — 0. For this kine- 
matical configuration, i.e. relative momentum in the Light-Front plane, the 
relativistic wave function ^f(q,9 = 90°) at high momentum is even found to 
be bigger than non relativistic one. To get rid of this dependence, we compared 
\^(q, 6)\ 2 integrated over the degree of freedom, both for relativistic and non 
relativistic solutions. The resulting functions, displayed in figure 4, measure 
the effective relativistic effects in the wavef unctions. At q = they remain at 
the level of 5%, once the energy is readjusted. In the high momentum region 
the relativistic solution is - as expected - smaller than the non relativistic one, 
but their differences reach a factor three at q — 2, and this for a moderate 
value of the coupling constant a = 0.5. 



mq)l 2 
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Fig. 4. Squared modulus of the wave function integrated over 9. The LFD solu- 
tion (solid line) and the reajusted non relativistic one (dot-dashed line) have same 
binding energy. The small top right graph is a zoom of the high momentum region. 

In the case /i ^ there exists a critical a below which there is no bound 
solution. Figure 5 represents the binding energy B as a function of the cou- 
pling constant a for different values of /i. They are compared with those pro- 
vided by BS equation in the same ladder approximation, whose kernel incor- 
porates higher order intermediate states. We have solved this equation using 
the method described in [17]. A first remark in comparing both approaches is 
that their results are seen to be close to each other. This fact is far from being 
obvious - specially for large values of coupling constant - due to the differ- 
ences in their ladder kernel. A quantitative estimation of their spread can be 
given by looking into an horizontal cut of figure 5, i.e. calculating the relative 
difference in the coupling constant (q;l_fd _ a Bs)/ a LFD for a fixed value of 
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Fig. 5. Binding energy as a function of a for Fig. 6. Differences in the coupling 
different values of \x in LFD (solid) and BS constant as a function of \i for 
(dashed) approaches fixed values B. 

the binding energy. The results, displayed in figure 6 for B = 1.0,0.1,0.01, 
show that relative differences (i) are decreasing functions of \i for all values 
of B (ii) increase with B but are limited to 10% for the strong binding case 
B = m which involves values of a > 5. This indicates the relatively weak 
importance of including higher Fock components in the ladder kernel even for 
strong couplings, as already discussed in [13]. 



It is interesting to study the weak binding limit of both relativistic approaches 
and compare them with the non relativistic calculations in the case fi ^ 0. The 
results are given in figure 7a for fx—1. They show on one hand that LFD and 
BS (solid lines) converge towards very close, though slightly different, values of 
the coupling constant (^ ~ 0.01). On the other hand one can see, contrary to 
the fx — case in figure 2, a dramatic departure of both relativistic approaches 
from a non relativistic theory (dot-dashed line), even for negligibles values of 
binding energy. The differences increase with fx as shown in figure 7b in which 
LFD and BS results are not distinguished. The origin of this departure lies 
in the fact that the integral term in equation (1) is dominated by the region 
q' ~ fx, even for very small values of B, and for the case fx ~ m the ^ terms - 
which make the difference between the non relativistic and relativistic kernels 
- are not longer negligible. We conclude from that to the non adequacy of a 
non relativistic treatment in describing systems interacting via massive fields, 
what is the case of all the strong interaction physics when not described via 
gluon exchange. 

Some approximations of equation (1) have been studied in order to disen- 
tangle the different contributions to the relativistic energies B(a) (see figure 
8). Equation (1) is formally written K^/ = j-V^. We first consider the case 
of a non relativistic kernel V - i.e. a Yukawa potential - and e q = m, with 
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(a) (b) 

Fig. 7. Zero binding energy limit of LFD and BS equations (solid) compared with 
the non relativistic solutions (dot-dashed) for fi = 1 (a) and for different values of 
M (b). 

in curve a the non relativistic kinematics K = 4q 2 + 2mB and in curve a' 
the relativistic one K = A(q 2 + m 2 ) — M 2 . Curves b and b' are obtained in 
the same manner, but putting e q = \Jq 2 + m 2 . The last one corresponds to 
the LFD equation. The results in figure 8 show that the kinematical term K 
has a very small influence on B(a), whereas the contributions of e q and V 
to the total binding are both essential. We conclude from this study that the 
kinematical corrections alone, as they are performed e.g. in minimal relativity 
calculations, are not representative of relativistic effects. Even by including 
them in the kernel through e q the results obtained are wrong by a factor 2. 

In case of an energy dependent kernels the normalization condition (5) is 
only approximate. This energy dependence denotes coupling to higher Fock 
components and the correct normalization condition for the model considered 
reads + 7V"( 3 ) = 1 where iV® is the norm contribution from the three-body 
Fock component. Using (5) only the two-body part is included. One can show 
that the correction to the two-body normalization condition is given by: 

^ (3) = -7^/— — « r V,? , -n)^(? ) 9', n,M 2 Mq,q-n) (6) 

(Z7r) D J E q E q i OM Z 



This expression can be analytically integrated over two angles ip and ip' and 
we are let with a four dimensional integration. The three-body correction to 
the norm, i.e. the ratio — 7% — , as a function of the coupling constant is 
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Fig. 8. Relativistic contributions of different terms in equation (1) (see text) 

shown in figure 9 for the case \x = 0.15. We remark that this correction is not 
zero at the critical value a = 0.35 corresponding to the B = threshold. Its 
behaviour in the region of large coupling tends asymptotically towards a value 
non exceeding 30%. This is in contrast with the evolution of the parameter 
R introduced in the same figure to estimate the norm correction for a system 
with a given value of R. For deuteron, e.g., one has R ~ 1% and the expected 
normalization corrections are of the order of 4%. 




Fig. 9. Normalization correction due to 3-body Fock components and relativistic 
parameter R = << L{' as functions of a for the case /x = 0.15. 
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4 A scalar model for Deuteron 



A simple relativistic model for deuteron in the LFD is obtained by adding 
to the interaction kernel (2) a repulsive part exactly analogous except for the 
sign of the coupling constant. Even if this procedure is no longer based on 
field theory - a scalar exchange cannot produce a repulsive interaction - the 
potential obtained constitutes a LFD relativistic version of the Malfliet and 
Tjon NN potential [18] on the form: 



The non relativistic model corresponds to Q 2 = (q — q') 2 and for the 3 Si 
state the parameter set Vr—7.29, Va=3.18, ^=0.314 GeV, \i#=1\ia inserted 
in Schrodinger equation ensures a deuteron binding energy B = 2.23 MeV. 

By solving the LFD equation (1) with potential (7) we can estimate the mod- 
ification in the deuteron description due to a fully relativistic treatment. The 
first result concerns its binding energy which becomes B = 0.96 MeV. The 
inclusion of relativity produces thus a dramatic repulsive effect, already drawn 
in [19]. We emphasize that, as mentioned before in the case of Wick-Cutkosky 
model, the use of relativistic kinematics alone induces a very small change 
in the binding energy. The sizeable energy decrease is almost entirely due to 
the r.h.s. part of (1). To obtain a proper deuteron description in a relativistic 
frame it is necessary to readjust the parameters of the non relativistic model. 
A binding energy of 2.23 MeV is recovered with a repulsive coupling constant 
Ar = 6.60 MeV - all other parameters being unchanged - what represents 
an decrease of 10% with respect to its original value. Another possibility is 
to increase the attractive coupling constant up to Va = 3.37. The relativistic 
effects in deuteron wave function depend sensibly on the way the energy is 
readjusted as well as on the relative angle 9 between n and the momentum 
q. For instance when modyfing Va and for the value 9 = 0, the zero of the 
relativistic wave function is shifted by m 0.1 GeV/c towards smaller values of 
q and the differences in the momentum region of q = 1.5 GeV/c are 50% in 
amplitude. 



5 Conclusion 



We have obtained the solutions for a scalar model in the Light Front Dynam- 
ics framework and in the ladder approximation. The results presented here 
concern the S-wave bound states. 
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We have found that the inclusion of relativity has a dramatic repulsive effect 
on binding energies even for systems with very small values. The effect is 
specially relevant when using a scalar model for deuteron: its binding energy 
is shifted from 2.23 MeV down to 0.96 MeV. This can be corrected by de- 
creasing of 10% the repulsive coupling constant, what indicates the difficulty 
in determining beyond this accuracy the value of strong coupling constants 
within a non relativistic framework. 

Light-Front wave functions strongly differ from their non relativistic counter- 
parts if they are calculated using the same values of the coupling constant. 
Once the interaction parameters are readjusted to get the same binding energy 
both solutions become closer but their differences are still sizeable. 

The relativistic effets are shown to be induced mainly by the relativistic terms 
of the kernel. The relativistic kinematics has only a small influence on the 
binding energy; furthermore, its effect is attractive whereas the total relativis- 
tic one is strongly repulsive. 

The normalization corrections due to the three-body Fock components increase 
rapidly for small values of the coupling constant and saturate at ~ 25% in the 
ultra relativistic region. They have been estimated to 4% in the deuteron. 

The LFD results are very close to those provided by Bethe-Salpeter equation 
for a wide range of coupling constant despite the different physical input in 
their ladder kernel. However in the case of systems interacting via a massive 
exchanged field, they both strongly differ from the non relativistic solutions 
even in the zero binding limit. This leads to the conclusion that such systems 
cannot be properly described by using a non relativistic dynamics. 

The case of higher angular momentum states for scalar particles requires more 
formal developments and is presented in a forthcoming publication [20]. 
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